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Abstract 

We examine the response of closed-shell nuclei using a correlated interaction, derived with the Unitary Correlation 
Operator Method (UCOM) from the Argonne V18 potential, in second RPA (SRPA) calculations. The same correlated 
two-body interaction is used to derive the Hartree-Fock ground state and the SRPA equations. Our results show that 
the coupling of particle-hole states to higher-order configurations produces sizable effects compared with first-order 
RPA. A much improved description of the isovector dipole and isoscalar quadrupole resonances is obtained, thanks 
in part to the more fundamental treatment of the nucleon effective mass offered by SRPA. The present work suggests 
the prospect of describing giant resonance properties realistically and consistently within extended RPA theories. 
Self-consistency issues of the present SRPA method and residual three-body effects are pointed out. 
Keywords: Second RPA, Giant Resonances, Unitary Correlation Operator Method, realistic effective interactions 
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The problem of deriving global effective nucleon- 
nucleon interactions (NNI), based on realistic NNIs 
and appropriate for use in microscopic many-body 
theories of the nucleus, has been addressed recently 
in the framework of the Unitary Correlation Opera- 
tor Method (UCOM) [1]. One starts from a realistic 
NNI and treats explicitly the short-range correla- 
tions it induces, so that a softened, phase-shift equiv- 
alent NNI is obtained. Such correlated realistic NNIs 
offer the possibility to exploit microscopic theories, 
ranging from Hartree-Fock (HP) and many-body 
perturbation theory to, e.g., no-core shell model, in 
a consistent and systematic way, as has been demon- 
strated in a series of applications [2|3|4j . 

The mean-field picture of the nucleus is empiri- 
cally justified and, thanks to its simplicity, it is the 
most convenient starting point to the microscopic 
description of nuclear response throughout the nu- 
clear chart. Small- amplitude oscillations of the den- 
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sity induced by an external field, in particular giant 
resonances (GRs), can be described self-consistently, 
within the Random Phase Approximation (RPA) 
(or quasiparticle RPA - QRPA) based upon the HP 
ground state (or HF Bogolioubov - HFB - to in- 
clude pairing). One can go beyond first-order RPA 
and describe coUisional damping and spreading of 
GRs within, e.g., second RPA (SRPA), or equiva- 
lent methods; higher-order effects are related to non- 
trivial parts of the two-body density matrix and can 
influence the position, strength and fine structure 
of GRs, as well as those of less collective low- lying 
states - see, e.g., Refs. |5l6l7l8l9ll0lllll2fT^ and 
references therein. 

The input to such calculations is an effective 
nuclear Hamiltonian or a density functional. Phe- 
nomenological effective NNIs are fitted to sets of 
experimental data using mostly IIF(B) calculations 
and selected (Q)RPA results. Their range of ap- 
plicability is inevitably restricted to the selected 
observables and many-body methodJ^. 



^ Let us note that, since no effective interactions have been 
fitted to SRPA calculations, mainly for computational rea- 
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The UCOM offers an alternative path to such mi- 
croscopic calculations. Applications in HF and RPA 
were discussed in Refs. [314114] . In this work we 
will use the same UCOM interaction, derived from 
the Argonne V18 potential, in SRPA calculations 
of nuclear response. As we will see, SRPA and the 
UCOM make an interesting combination: Not only 
can SRPA accommodate more physics than first- 
order RPA, it also appears suitable for describing 
long-range correlations (LRC) which are excluded 
from the UCOM by construction. 

It is instructive to summarize the properties of the 
UCOM and the main conclusions of previous work. 
Short-range correlations, both central and tensor, 
are explicitly described within the UCOM by means 
of a unitary correlation operator. This can be used 
to perform a similarity transformation of the bare 
nuclear Hamiltonian (or any other operator of in- 
terest). The resulting transformed, or "correlated", 
Hamiltonian consists of a kinetic-energy term and 
an energy-independent two-body potential Vucom • 
The UCOM potential is phase-shift equivalent to the 
original, bare one and has been shown to have good 
convergence properties 2 3 . Omitted three-body ef- 
fects (correlations and interactions) are effectively 
taken into account to some extent by the parameter- 
ization of the correlators, while the task of describ- 
ing LRC is assigned to the model space. The only pa- 
rameters entering the formalism are in fact related to 
the optimal separation of state-independent short- 
range correlations and the minimization of three- 
body effects. They are fixed in the nucleus "'He f?. 

The VucOM was employed in HF calculations in 
Ref. [3]. The binding energies and charge radii are 
underestimated in HF and the level spacing of the 
single-particle states is too large. Second-order per- 
turbation theory (PT) constitutes a seemingly ade- 
quate extension of HF [3] . A very good description 
of nuclear binding energies was achieved within PT 
for nuclei from ''He to ^"*Pb, suggesting that the HF 
underbinding (about 4 MeV per nucleon) can be at- 
tributed to missing LRC. Charge radii are still un- 
derestimated within PT, implying that supplement- 
ing the two-body Vucom Hamiltonian with a three- 
body term (and readjusting it accordingly) to take 
into account missing effects may be necessary. 

The VucoM was subsequently employed in stan- 
dard, HF-based, self-consistent RPA calculations 



sons, consistency in the treatment of the ground and excited 
states is, more often than not, lost in practical applications 
beyond RPA. 



to study nuclear GRs [3]. The isoscalar (IS) giant 
monopole resonance (GMR), the isovector (IV) 
giant dipole resonance (GDR), and the IS giant 
quadrupole resonance (GQR) were examined. A 
reasonable agreement with the experimental GMR 
centroid energies was achieved for various closed- 
shell nuclei, but the energies of the GDR and the 
GQR were overestimated by several MeV. 

High GDR and GQR energies, as well as low 
single-particle level densities compared with exper- 
iment, are usually associated with a small nucleon 
effective mass, which in the case of the Vucom 
can be viewed as a result of missing residual LRC. 
Including explicit RPA correlations in the ground 
state was shown to produce small corrections [Tl] . 
Given that an extended model space is of great rel- 
evance when using the Vucom, it is important to 
examine whether coupling of the particle-hole (ph) 
excitations to higher-order configurations {2p2h 
and beyond), starting with SRPA, can produce 
significant corrections. 

We will use the SRPA as it was formulated in 
Ref. [ini in analogy to RPA. Excited states |A) of 
energy Ex = hujx with respect to the ground state 
|0) are considered as combinations of ph and 2p2h 
configurations. The corresponding creation opera- 
tors Q\, such that 

|A)-Qi|0), Qa|0) = 0, (1) 
are then written as 



'^Plhip2h2'^PlhiP2h2 '^pihip2h 
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where O^j^ creates a ph state and 0^^^,^, creates a 
2p2h state. We omit angular momentum coupling to 
keep the notation simple. The SRPA ground state, 
which formally is the vacuum of the annihilation 
operators Qx, is approximated by the HF ground 
state. The forward {X, X) and backward (F, y) 
amplitudes are the solutions of the SRPA equations 
in ph 2p2/i— space 
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where A and B are the usual RPA matrices, A12 de- 
scribes the couphng between ph and 2p2h states and 
A22 contains the 2p2h states and their interactions. 
If we neglect the coupling amongst those states, A22 
is diagonal and its elements are determined by the 
unperturbed 2p2h energies, 

[A22]pihip2h2,p[h'^p'2h',^ ~ 
^Pip'i^hih'^^pip'^^hih'^i^pi "I" Chi S^ij), (4) 

where are the HF single-particle energies. It is 
interesting to note that the SRPA problem of Eq. ^ 
can be reduced to an energy-dependent eigenvalue 
problem of the dimension of the RPA matrix (see, 
e.g., Ref. [6]). Therefore, it can be viewed as an RPA 
problem with an energy-dependent interaction. 

The dimension N of the SRPA matrix, Eq. ©, 
can be rather large. For the purposes of the present 
work we solve problems with N up to 10^, but larger 
spaces are to be expected for heavier nuclei and 
larger bases. Fortunately, the SRPA matrix is also 
sparse, especially, but not only, when the approxi- 
mation (HI) is employed. Thus it is possible to store 
all its non-zero matrix elements and then use a Lanc- 
zos procedure to obtain only the spectrum section 
of interest, i.e., a couple of hundred eigenstates and 
eigenvalues at the lower end of the spectrum. 

By setting the coupling matrices A12 and A22 and 
the 2p2h amplitudes X, y equal to zero in Eq. dS]), 
the usual RPA problem is retrieved. If in addition we 
neglect the ph residual interaction (i.e., Bph^pih' = 
and Aph,p>h> = (cp - eh)6pp'dhh'), a trivial, unper- 
turbed problem is obtained, where the eigenstates 
I A) are simply the ph configurations at the HF level 
and the Y amplitudes vanish. In all cases we may 
define the strength distribution of a \ph~^) configu- 
ration as the quantity 

Spf,{Ex) ^ \X^,,\' ~ \Y^\\\ (5) 

Note that in the unperturbed case the centroid of 
Sph{E) is identical to the ph energy Cp — e/j and its 
width is zero. 

As usual, the quantities of interest are transition 
strength distributions Rf{E) of single-particle op- 
erators F'f = Y,ij fij^laj, 

Rf{E) = Y,\{MF^O)\^S{E-Ex) (6) 

A 

= Y,Bf{Ex)S{E-E^), (7) 

A 

and their energy moments 
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Fig. 1. IV dipole response of ^'^O in a single-particle basis 
with (2n + <?)max = 12 and £max = 8. Pale (cyan) dash- 
dotted lines: full SRPA solution; dashed (black) lines: using 
approximation l|4]l; full (red) lines: RPA. 

mfc = ^ii;^i?;^(i?A), (8) 
A 

determined by the amplitudes X and Y through 
(A|Ft|0) - Y^^fphX^^* + hpYX]- (9) 

ph 

We will consider IS and IV transitions of definite 
spin and parity , described by standard single- 
particle transition operators [4'. 

The total strenght mo and the first moment of 
the strength distribution mi are the same in SRPA 
as in RPA ^TE\. However, when based on the HF 
ground state, the SRPA is not fully self-consistent 
and symmetry-conserving, contrary to the HF-based 
RPA, as it misses a class of second-order effects re- 
lated to ground-state correlations |17I18| . The miss- 
ing effects may be important, especially for the less 
collective low-lying states. In principle, it is possi- 
ble to combine the SRPA with a correlated ground 
state |7I17I19| . or employ a self-consistent Green's 
function method '20] , for a more complete theoreti- 
cal treatment of nuclear excitations, but that is be- 
yond the purposes of the present work. 

We have used the same Vucom as in Refs. [21314] 
and a harmonic-oscillator single-particle basis 
(length parameter b — 1.7fm) and we have exam- 
ined the IS monopole (ISM), IV dipole (IVD) and 
IS quadrupole (ISQ) response mainly of the nuclei 
^^O and '*'^Ca in SRPA. The diagonal approxima- 
tion, Eq. ^ is used. It has been verified, though, 
that inclusion of the 2p2h couplings does not intro- 
duce significant corrections. An example is shown 
in Fig. [TJ Note that those couplings constitute 
higher-order effects and their smallness suggests 
that corrections beyond second order are not large. 

In what follows, single-particle states with radial 
quantum number up to rtmax — 6 and orbital angular 
momentum up to ^max = 6 have been included. The 



3 




1 


1 1 1 1 

I RPA - 




SRPA - 




■\ '°Ca - 




1 .■; ISM _ 
\ * * 











40, 



Ca 
ISQ 



•I 1 - 



10 



20 30 40 
E [MeV] 



50 





1 «l 


1 1 1 

experiment o _ 




° l\ '■ 








^°Ca - 















60 



Fig. 2. The IS monopole (top), IV dipole (middle) and IS quadrupole (bottom) strength distributions for the nuclei ^''O (left) 
and '^'^Ca (right) within RPA (blue dashed lines) and SRPA (red full lines), compared with experiment (points, arrows). The 
calculated distributions have been folded with a Lorenzian with a width of 2 MeV. The experimental centroids mi /mo of the 
ISM and the ISQ GRs were taken from Refs. [2l] (1*^0) and [22] (4°Ca). Photoabsorption cross sections were found in Refs. 
| 23I24| i^^O) and [25] (*°Ca)(data available in [2g]) and experimental IVD strength distributions were evaluated from those. 



convergence of the GR sum rules mo and mi and 
centroids is rather good for the present basis (within 
about 1-2 MeV for the centroids, being worst for the 
IS GMR). 

The spurious state related to the CM motion will 
generally not be exactly seperated from the physical 
spectrum, when SRPA is based on the HF ground 
state [18| . In order to quantify this problem, we have 
examined the IS dipole response. We found that spu- 
rious states appear at about 5-8 MeV. Using a transi- 
tion operator of the usual radial form (oc r'^ — |(r^)r) 
and its uncorrected form (oc r'^), we found that the 
spectrum in the GR region is not strongly affected 
by the choice of operator and can be considered un- 
contaminated. Further technical and numerical de- 
tails regarding our SRPA implementation and con- 
sistency tests are reserved for a more extended fu- 
ture publication. 

In Fig. m we show the ISM, IVD and ISQ strength 
distributions for ^^O and ''"Ca. The calculated 
spectra (RPA and SRPA) have been folded with a 
Lorenzian with a width of 2 MeV, for presentation 



purposes. In all cases, we observe that the SRPA 
centroid energies are much lower than the RPA 
ones. The reason for the large difference between 
the RPA and SRPA results - even for such collec- 
tive ph excitations like the GRs considered here - 
is, to a large extent, the coupling of ph states with 
virtual phonons, implicitly taken into account in 
SRPA. The inclusion of second-order configurations 
within SRPA effectively dresses the underlying 
HF single-particle states with self-energy inser- 
tions [5|6|8j and brings them closer to each other 
energetically, thereby lowering the ph energies. As 
an illustration of the effect, we show in Fig. |3|how 
the strength of the dipole configurations corre- 
sponding to a neutron 0^3/2 hole and a neutron 
np3/2 particle {n = 1,2,...) is distributed in HF, 
RPA, and SRPA. In HF all the strength of each ph 
configuration, defined in Eq. ([5]), is concentrated 
in one peak, positioned at energy equal to Cp — eh- 
In RPA the strength of each configuration appears 
slightly shifted and fragmented. In SRPA the shift 
and fragmentation are much more pronounced. The 
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Fig. 3. Fragmentation and shift o{ ph states. Bars show how 
the spectroscopic strength Sph{E), Eq. ((Sj, of the ph config- 
urations \{i>Ps/2){'^0'i3/2)~^^ 1~) (contributing to the dipole 
strength) is distributed in ''"Ca, within HF (bottom), RPA 
(middle) and SRPA (top). Thicker, pale (cyan) bars denote 
the distribution of \(i^^P3/2)('-'0d^/2)~^ 't 1~) (one shell) only. 
Continuous (blue) lines, corresponding to the y— axis scale 
on the right: IVD strength distribution, shown for reference. 

shift is related to the real part of the acquired self 
energy and the fragmentation to the imaginary part 
and neither can be ignored when using completely 
"undressed" (with respect to LRC) HF states like 
the ones produced by the Vucom- In this scheme 
the HF single-particle energies (in a similar man- 
ner as HF binding energies and radii) are viewed 
as auxiliary model quantities, not to be directly 
compared with experimenif^. Double counting of 
second-order effects is thus avoided. 

In Fig. [2] we also observe that, in some cases, the 
resonance width seems to be smaller in SRPA than 
in RPA, contrary to what one might expect. The 
reason appears to be again an overal compression of 
the underlying ph spectrum that competes with the 



^ A related discussion in the context of the Extended Theory 
of Finite Fermi Systems and phenomenological effective NNIs 
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Fig. 4. Fragmentation in SRPA. Both panels show the IVD 
transition strength distribution of *'^Ca in RPA (black bold 
bars) and SRPA (pale cyan bars). In the lower panel the 
y— scale is different, in order for the large amount of weaker 
SRPA states to become visible. 
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Fig. 5. IVD transition strength distribution of ^'^Zr in RPA 
and SRPA compared with (■y,Xn) data from |28I26I . 

increased fragmentation. Fragmentation does occur 
in SRPA, as is illustrated in Fig. [H 

Let us now discuss our results in comparison with 
experiment. In the middle panels of Fig. [2] our IVD 
strength distributions are shown along with those 
extracted from experimental data (there has been 
no ad hoc renormalization imposed). We observe 
that the IV GDR is more realistically reproduced by 
SRPA than by RPA. Its centroid energy is somewhat 
underestimated. These trends seem to persist in the 
heavier nucleus ^"Zr, as shown in Fig.[51 In the lower 
panels of Fig. [2] we show the ISQ strength distribu- 
tions. The RPA and SRPA results are shown and the 
experimental centroids of the IS GQR are indicated. 
The SRPA results agree very well with experiment, 
suggesting that the coupling to higher-order config- 
urations restores a realistic nucleon effective mass 
to a large extent. For the distribution widths, which 
we evaluate from the energy moments in the reso- 
nance region as cr = \/m2/mo — (mi/mo)^, we find 
3.8 MeV (2.5 MeV) for ^^O C^Ca), smaller than, 
but comparable with, the experimental widths of 



5 



5.1 MeV (2.9 MeV) [HMj- SRPA is a suitable the- 
ory for describing the spreading width of resonances, 
but one should keep in mind the important role of 
the continuum in light and medium-heavy nuclei, 
which is not properly described by our method and 
model space. 

In the upper panels of Fig. [2] the ISM strength 
distributions are shown. The energies of the IS 
GMR are underestimated within SRPA. This may 
be another indication that there are missing three- 
body effects. Normally, residual three-body correc- 
tions should affect the IS GMR most of all, since 
it is a compression mode. They should affect less 
strongly the IV GDR, where the nuclear interior 
plays a lesser role, and less the IS GQR, which is 
a surface mode. These physical arguments could 
serve as a guide for the construction of an appro- 
priate effective three-body term to supplement the 
two-body VucoM • In any case, the three-body term 
shall depend on the bare interaction used, as well as 
the correlation operators applied to it [2] . Another 
source of problems could of course be the inher- 
ent inconsistencies of the present SRPA formalism, 
which will be investigated in future work. 

In summary, we have used a correlated interac- 
tion, derived within the UCOM from the Argonne 
VI 8 potential, in SRPA calculations of nuclear re- 
sponse. Our results for the nuclei ^^O and ^°Ca 
show that the coupling to higher-order configura- 
tions produces sizable effects, compared with first- 
order RPA. An improved description of the IVD and 
ISQ resonances is obtained. Through a more funda- 
mental treatment of the nucleon effective mass, the 
UCOM-based SRPA method seems to enable a si- 
multaneous description of IVD and ISQ GRs. Our 
correlated interaction underestimates the energy of 
the ISM GR, though, pointing to missing three-body 
effects. Efforts to construct appropriate three-body 
potentials to complement our two-body correlated 
potentials are under way. 

The present work suggests the prospect of de- 
scribing GR centroids and structure realistically 
and consistently, within extended RPA theories 
like SRPA. The crucial points in this context are 
that the Vucom interaction does not parameterize 
LRC, which are instead described by the SRPA, 
and of course the good convergence properties of 
the Vucom, which should render the SRPA (or 
equivalent) model space flexible enough to describe 
residual correlations. More systematic calculations 
are planned for the immediate future, in order 
to assess and explore this possibility. These shall 



include low-lying collective states. Resolving the 
self-consistency issues of the present SRPA method 
is planned as well. 
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